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Abstract 

We study the issues of scaling and universality in spectral and trans- 
port properties of the infinite dimensional particle-hole symmetric (half- 
filled) Hubbard model within dynamical mean field theory. One of the 
simplest and extensively used impurity solvers, namely the iterated per- 
turbation theory approach is reformulated to avoid problems such as ana- 
lytic continuation of Matsubara frequency quantities or calculating multi- 
dimensional integrals, while taking full account of the very sharp struc- 
tures in the Green's functions that arise close to the Mott transitions 
and in the Mott insulator regime. We demonstrate its viability for the 
half-filled Hubbard model. Previous known results are reproduced within 
the present approach. The universal behavior of the spectral functions 
in the Fermi liquid regime is emphasized, and adiabatic continuity to the 
non-interacting limit is demonstrated. The dc resistivity in the metallic 
regime is known to be a non-monotonic function of temperature with a 
'coherence peak'. This feature is shown to be a universal feature occurring 
at a temperature roughly equal to the low energy scale of the system. A 
comparison to pressure dependent dc resistivity experiments on Selenium 
doped N1S2 yields qualitatively good agreement. Resistivity hysteresis 
across the Mott transition is shown to be described qualitatively within 
the present framework. A direct comparison of the thermal hysteresis ob- 
served in V2O3 with our theoretical results yields a value of the hopping 
integral, which we find to be in the range estimated through first-principle 
methods. Finally, a systematic study of optical conductivity is carried out 
and the changes in absorption as a result of varying interaction strength 
and temperature are identified. 

Keywords: Hubbard model; dynamical mean field theory; iterated per- 
turbation theory; optical conductivity; hysteresis 
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1 Introduction 



The Hubbard model (HM) [TJ [U [3j 0] has become one of the central paradigms 
in theories for strongly correlated electron systems. For systems ranging from 
ultracold gases [5] , transition metal oxides [5J [7] , high temperature supercon- 
ductors [SI H] , and even graphene [TUJ E] , the HM has been employed in various 
guises. Thus, the importance of developing an understanding of this model 
cannot be overstated. The theoretical techniques that have been employed 
in an attempt to solve this model are diverse, e.g. Hatree-Fock approxima- 
tion [12l[T3l[T4], Hubbard's Green's function approximations [HI [16], variational 
wave function method [4[ [17], slave Boson approach [18], etc. However, to date 
the general HM remains unsolved except for the one-dimensional case [19] . A lot 
of progress has been made in recent years in this direction through the dynami- 
cal mean field theory (DMFT) [20, 21], where the lattice model is mapped onto 
an effective single impurity model with a self-consistently determined hybridiza- 
tion. Within the DMFT context, a greater interest is being devoted towards 
the inifinite dimensional HM since the mapping becomes exact in the limit and 
simplifies many major issues of many-body approaches. Although the lattice 
problem reduces to that of an impurity problem embedded in a non-interacting 
bath, the problem remains amenable to a host of techniques such as non-crossing 
approximation (NCA) [22], quantum Monte Carlo(QMC) [HO [21], numerical 
renormalization group (NRG) [24] , exact diagonalization (ED) [25] , iterated per- 
turbation theory (IPT) gFJ [23 [13 [53] , local moment approach (LMA) [2"§I|gO"] . 
etc, that have been developed to solve the impurity problem. IPT is a partic- 
ularly simple diagrammatic perturbation theory based approach that has been 
extensively used to study not only the impurity problem, but also a host of 
other lattice based models such as the HM and the periodic Anderson model 
(PAM) [21]. It is known to yield qualitatively good results, such as the Mott 
metal-insulator transition, heavy fermion behavior etc. However, since it is 
perturbative by construction, IPT does not do as well in quantitative terms. 
The half-filled HM, has likewise, been a very intensively studied model, as a 
description of the Mott MIT, and IPT has been very successful in elucidating 
quite a few conceptual issues (27] [28]. Nevertheless, we observe a few important 
deficiencies in the present implementations of IPT, specifically for the half-filled 
HM. IPT, as implemented generally, calculates Matsubara frequency Green's 
functions and self-energies. The procedure of analytic continuation, required to 
obtain real frequency quantities, is known to be a numerically difficult prob- 
lem, which is conventionally treated by Pade approximation or the maximum 
entropy method |32[ I21j . The direct real frequency implementation takes two 
routes - in the first, one has to solve multi-dimensional integrals, which is com- 
putationally expensive; the other route employs Fast Fourier transforms, which 
has the disadvantage that the frequency grid is necessarily uniform, and thus 
cannot sample the very sharp structures arising in the Green's functions in the 
parameter region proximal to the Mott transition as well as in the Mott insu- 
lator regime. This leads to problems in convergence of the DMFT iterations. 
In a very recent study, an ad-hoc way of dealing with the sharp structures was 



2 



implemented through the addition of a small constant imaginary part added to 
the real frequency [33]. This method broadens the poles of the Green's functions 
and lets one handle the sharp structures numerically. However, it has the dis- 
advantage that it leads to spurious spectral features, e.g. the insulating phase 
does not exhibit a clear gap. 

In this work, we introduce an alternative real frequency implementation that 
overcomes the above limitations, and hence, lets one study the properties of 
interest such as spectral functions etc directly on the real axis with an adaptive 
non-uniform grid. In order to solve the DMFT self-consistency equations we 
start with an initial guess for the self-energy (typically the Hatree-Fock) whereas 
the host Green's function is guessed initially in common practice. We benchmark 
our results against known IPT results calculated through analytic continuation 
of the Matsubara frequency Green's functions. We find excellent agreement. 
Issues of universality and scaling of spectral functions and transport quantities 
are emphasized. The paper is organized as follows: The next section details 
the model and the reformulation of IPT. The numerical implementation details 
are also discussed. Section 3 comprises our theoretical results, discussion and 
comparison to experiments. We present our conclusions in section 4. 



2 Methods and Formalism 

The single-band HM in standard notation is given by 

H=- ^ U j {cl r7 c ja +h.c.) + e d J2 c L c ta + u Yl n ^ n ^ W 

<ij>,cr i<7 i 

where Uj is the amplitude of hopping of electrons from site-z to site-j , e d is the 
electron's orbital energy and U is the on-site Coulomb repulsion. The operator 
Ci„ {c\ a ) annihilates (creates) an electron of spin a at site i. In the limit of 
infinite dimensions (d — > oo), the hopping parameter is scaled as t — t^/yfz 
where z is the coordination number of the lattice. In this work, we choose to 
work with the hypercubic lattice (HCL) for which the non-interacting (U = 0) 
density of states is an unbounded gaussian, i.e. 

Do (e) = -L-cx P (-eytl) (2) 

The major simplification within DMFT, which is exact in d — oo, is that the self- 
energy [20j and vertex function |34j become momentum independent or spatially 
local. The local, retarded Green's function in the paramagnetic case is given by 

where ui + = uj + irj, rj — >• + and S(o;) is the real-frequency self-energy. The k 
sum may be transformed to a density of states integral and thus may be written 

as 

G{lj) =n[y(u)] (4) 
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where H[z] = J de Do(e) / (z — e) is the Hilbert transform of z with respect to 
Do(e) and 7(w) = uj + — td — Since a lattice model can be mapped onto 

an impurity model within DMFT with the self-consistency condition that the 
impurity self-energy is the same as the lattice self-energy, one can find the host 
or medium Green's function G(uS) for the impurity through the Dyson equation 



(5) 



Solution of the impurity model in terms of G(u)) would then yield a new 
which when put back in Eq. Q gives G(uS). Thus, given an impurity solver 
technique, one can self-consistently solve for the Green's functions, and hence 
the self-energy. Since the vertex function is also momentum independent within 
DMFT, the calculation of conductivity through the current-current correlation 
function involves only single-particle quantities. Thus the optical conductivity 
may be calculated by using the following expression obtained from the Kubo 
formula 1351. 
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where <jq = Aite 2 t 2 a 2 n/K for lattice constant a, electronic charge e, and electron 
density n. As uj — > we obtain the DC conductivity : 
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where D{uj) = — jlmG(w) is the spectral function of the retarded Green's func- 
tion G(oj) and n F {uj) = (e^ + l)" 1 is the Fermi function (/9 = \/{k B T)). 



As mentioned earlier, the full solution of DMFT problem requires an im- 
purity solver technique. In this paper, we have chosen IPT for solving the 
self-consistent impurity problem. IPT does have a few drawbacks. It is based 
on a simple truncation of the diagrammatic perturbation theory in U to second 
order around the Hartree mean field solution. Thus it needs to be benchmarked 
for every new problem. For the same reason, it is unable to capture exponen- 
tially small energy scales. Nevertheless, the reason for this choice is that IPT is 
technically one of the simplest approaches that can capture the Mott transition 
physics in a qualitatively correct manner. IPT has been extensively employed 
by various groups to study the HM [571 [551 HI] and PAM [30 [37J [351 [35] . How- 
ever, the previous implementations suffer from the problems mentioned in the 
introduction section. 

We now describe IPT and our implementation of it briefly. The ansatz for 
the dynamical part of the self-energy (apart from the static Hartree term) within 
IPT is just the second order term of the perturbative expansion in U about the 
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Hartree limit, i.e. 

U 2 



S 2 (w) = lim + X! ^ ^ w + *^m)(? («^ P + w m )£ (iwp) (8) 



where iuj and it/ denote odd and even Matsubara frequencies respectively. Us- 

/oo j j (uj r *) 
diu' , (where Dg(uj) = 
-oo lUJ n -Uj' 

—\vaQ{u3)/'K) and carrying out the Matsubara sums along with the trivial ana- 
lytic continuation iu> n — > lj + , we get the following expressions for the imaginary 
part of the self-energy on the real frequency axis: 



De(uj) = ImE(w) 

7T 



U 2 / dw'D g (-u>') 



n F (u)')x(w + w') + n F (— uj')x{— u - w') 



0) 



where 



X(uj) = / dw'Dg(ui + ui')D g (u}')n F (-U)-U>')n F (u}') (10) 

J — OO 

and we get the real part by using Kramers-Kronig transformation: 

ReEM^f^^M (11) 

Eq. ([3]) , ((4j) , and ([5]) along with Eq. ([9l)- (fTT]) constitute the necessary ingredients 
for the solution of the half-filled Hubbard model within DMFT. We remark here 
that the integrals above are all one-dimensional and hence do not present any 
excessive computational expense. 

Numerical implementation of the above equations is straightforward in al- 
most the whole metallic regime. However, the metallic regime in proximity to 
the Mott transition and the Mott insulating regime is trickier. We will illustrate 
this point and its resolution for the Mott insulator here and for the correlated 
metal in the appendix. The retarded host Green's function may be separated 
into a singular and a regular part (C? reg ). 

g(u)=Y — +r g H, 77 ^0+, a>0 (12) 

i 

The sum in the above equation is over the poles uii or the singularities of the 
G(oj) and aj's are the corresponding weights or residues. 

It is easily seen that if 9(u>) has to satisfy the self-consistent equations of 
DMFT [([5|-(j5J).([5]l-([TT]l]. the residue a< of the singularities must satisfy self- 
consistent equations. For example, in the Mott insulator case where there is 
just one pole with residue a at the Fermi level (this follows from the singular 
behavior of the self-energy and host Green's function in the atomic limit (t* = 
0) [U EQlOn]), we get a cubic equation (at T = 0), given by 

«"' = 1 + IS < 13 > 

U or 
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where M 2 is the second moment of D (uj) [M 2 = duj uj 2 D {lo)]. Solutions 
of this equation by Cardano's method [41] shows that out of three roots, the 
physically reasonable root (such that a — > 1 as U — > oo) exists only for U > 
U c i ~ 3.67i* for the HCL. For a general U, we solve the above cubic equation 
numerically to get a. In the strongly correlated metallic regime, poles can 
arise in the host Green's function Q(io) even though the spectral density of 
the interacting Green's function shows a FL behavior within a width of low 
energy scale (described as u>l later). Two poles occur symmetrically about the 
Fermi level at zLluq proportional to the square-root quasiparticle weight Z (see 
Eq. (|24l0 . The Mott transition from the metallic regime to the insulating regime 
occurs through the collapse of these two poles at the Fermi level into one single 
pole characteristic of the insulating regime (as described above). The poles of 
the G(oj) lead to divergences at ±3wo in the (through Eq. flSJ). Using 

such a pole structure of the Green's functions and the self-energy, the critical 
U at which the metal transforms into the insulator, i.e U C 2 may be estimated 
as 4.77£* (see Appendix 2]) for the HCL. We also found the same for the Bethe 
lattice {U C 2 = 3.28t*) which agrees with the numerical estimation in Ref. . 

At finite temperatures, the singularity at the Fermi level must have a finite 
width. However, this width could be exponentially small in practice, and it is 
next to impossible to capture such sharp resonance numerically. So we utilize 
the spectral weight sum rule to compute the weight of the singularity, i.e. a. 
The presence of this singularity is numerically detected by a significant deviation 
(in practice 2%) of the integrated spectral weight of Dg from unity. The sharp 
resonance at the Fermi level in Dg is then numerically cut off to get T> r g S and 
the weight a is obtained by a — 1 — duTfg 5 . With the above separation 
of Dg into regular and singular parts, the self-energy expression reduces to 

U 2 a 3 

E( W ) - + __ (14) 

where E reg is obtained through the Kramers-Kronig transform of -D^ s which is 
given by 



D^{u) = U 



doj'D^^M^x^ + u:') + ^-n F {-uj)D^(u) + | ^(ui) 

dw' D g es (-u')n F (-u') X ies (-" - + ^ F H^ cg M + f X rog (-^) 

(15) 

and 

/oo 
dJ D Q c% (J)D g ^(uj+J)n F (J)n F (-u-J)+uD Q c% (uj)n F (-u) 
-OO 

(16) 

Now we proceed to discuss our results. 
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3 Results and discussions 



We have implemented the above mentioned formalism and have computed real 
frequency spectral functions, dc and optical conductivities. 



3.1 Spectral and transport properties 




Figure 1: (Color online) Spectral functions for various U showing quasiparti- 
cle resonances at Fermi level (left). A simple a;- axis re-scaling by Zt* leads 
to a collapse of all of these onto the non-interacting spectral function in the 
neighbourhood of the Fermi level thus signifying adiabatic continuity and FL 
behavior. 

Since our implementation is new, we would like to benchmark it by com- 
paring our results with other implementations. We begin with the local density 
of states (DoS) as given by D{uS) = — ImG(w)/7r. In Fig. [TJ we show the D(uj) 
computed at temperature T = for various U in the metallic regime. The left 
panel shows the spectral function on 'non-universal' scales, i.e. vs. w/t*, where 
the Hubbard bands are seen to form with increasing U, while the Abrikosov- 
Suhl resonance at the Fermi level gets narrower. The zero frequency value of the 
spectral function is also seen to be pinned at the non-interacting (U = 0) value. 
Such behavior finds a natural explanation in Fermi liquid (FL) behavior. A sim- 
ple linear expansion of the self-energy as = S(0) + (l — ^) oj + 0{uj 2 ) when 
used in Eq. Q gives G(w) = Gu=o{u/lul) where ujl — Zt* is the low energy 
Fermi liquid scale. Since the quasiparticle weight Z decreases with increasing 
U, thus signifying an increase in effective mass (m* = m/Z), the FWHM for 
a hypercubic lattice that is given by A = 2ojiy / ln(2) would decrease with in- 
creasing U . Another inference from such low frequency scaling behavior is the 
universality of the low frequency part of the spectral function. On the right 
panel of Fig. [TJ the same spectra as the left panel are plotted as a function of 
uj/ojl, and they are all seen to collapse onto the non-interacting limit spectra 
in the neighbourhood of the Fermi level. All of the above behavior is of course 
well known and well understood. We nevertheless emphasize that the extent 
of the FL regime is very small in strong coupling, because as the right panel 
shows, the scaling collapse is valid for uj < ljl, where ojl is expected to decrease 
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exponentially with increasing U. The approximation of IPT does not capture 
such exponential decrease of u>l, instead predicting an algebraic decrease. Nev- 
ertheless, the qualitative behavior from IPT indeed corroborates with other 
methods such as QMC. IPT has been shown to capture the first order Mott 
metal-insulator transition qualitatively and in reasonable agreement with other 
techniques such as NRG [21] . In the T — U plane, the metallic and insulating 
solutions are known to coexist in a certain region bounded by spinodal lines. 
The T = bounds are denoted by U c i and U C 2- This first order coexistence 
region may be simply found by computing the temperature or U dependence 
of the Fermi level density of states (-D(O) = — ImG(w)/7r) for various values of 
U/t* or T/t*. The resulting phase diagram is shown below in Fig. [2] It agrees 
well with those reported previously (2T] [28] . 
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Figure 2: The U — T 'phase-diagram' as calculated within IPT is shown here. 
The solid lines mark the spinodals of the coexistence region where the metallic 
and the insulating solutions coexist. 

The metallic region is a strongly renormalized FL, and thus the properties 
close to the Fermi level must be governed by a single low energy scale u>l- For 
universality to hold in the strong coupling region, the spectral function must 
have the following form: 

DfaT)=f(— ■— ) (17) 

i.e. it must be a function purely of T = T/cul and Q> = lo/lul [31j . In Fig. [3] 
we show the spectra in strong coupling for fixed T = 0.2 and increasing U/t*. 
We see that a scaling collapse does not occur implying that the above universal 
form does not describe the finite temperature IPT results. This non-universal 
behavior is an artefact of the specific iterated perturbation theory ansatz for 
the self-energy which is known to yield a non-universal form for its imaginary 
part [42], namely ImSi PT (w = 0) cx U 2 T 2 /tl- 
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Figure 3: Temperature scaling violated, as a function of temperature for U ft* — 
2.6 (dotted), 2.8 (dashed), 3.0 (dotted and dashed), 3.2 (doubly dotted and 
singly dashed), 3.4 (full line). 

The dc resistivity is computed using (|7|). Again, the pathologies in the 
imaginary part of the self-energy are reflected in the low temperature FL region 
in the following way. Although the resistivity does have a T 2 form, the expected 
universal scaling form of (T/u>l) 2 is not obtained [42]. However, surprisingly, the 
'coherence peak' position, which represents a crossover between low temperature 
coherent behavior to high temperature incoherent behavior does seem to be a 
universal feature in strong coupling as seen in Fig. U occurring at T = 0.6. 
The crosses represent the peak position, and as is seen in the right panel, these 
crosses line up at a single T. Thus the position of the coherence peak may be 
used to infer the low energy scale in a real material. In the inset of the left 
panel, we show the experimentally measured [43l|44] resistivity of NiS2-a;Se x as 
a function of pressure in Kbar (indicated as numbers). The resistivity for the 
lowest pressure rises dramatically with increasing T, before reaching a coherence 
maximum, and then decreasing slowly for higher temperatures. With increasing 
pressures, the initial rise becomes more gradual, and the coherence peak shifts 
to higher temperatures. An increase in pressure leads to a decrease in lattice 
spacing, thus increasing t*, the hopping parameter, while the local Coulomb 
repulsion U remains unaffected. Thus increasing pressure can be interpreted as 
a decrease in the U /t* ratio. A comparison of the inset with the main figure of 
the left panel clearly indicates qualitative agreement. The initial rise of p(T) 
with T is much sharper in experiment than in theory, but the rest of the features, 
including a shift of the coherence peak to higher temperatures with increasing 
pressure, are indeed observed. We emphasize here that the agreement is only 
qualitative and as such, no attempt is made to obtain quantitative agreement. 

We now study the resistivity hysteresis as obtained within IPT approxima- 
tion. 
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Figure 4: (Color online) Left panel: Theoretically computed resistivity as a 
function of temperature for various U/t* values (indicated as numbers). The 
crosses indicate the position of the coherence peak. Left panel inset: Experi- 
mentally measured [131 HI] resistivity for NiS2_ a; Se a ; as a function of pressure 
(indicated as numbers). Right panel: The same resistivities as in the left panel 
with the temperature rescaled by the low energy coherence scale u>l = Zt* 
showing that the coherence peak is indeed a universal feature of the strongly 
correlated metallic regime. 



3.2 Thermal hysteresis 

Fig. [5] shows the thermal hysteresis in resistivity obtained through heating and 
cooling cycles for fixed U's (i* = 1) in the coexistence region. The area enclosed 
by the hysteresis loop decreases as U — > U c . A full cycle hysteresis is observed 
only in the region f7 c (~ 3.5) < U < C/ C i(~ 3.7). For U < U c , the metal- 
insulator transition is continuous and hence of second order while for U > U c i 
the transition is discontinuous but hysteresis is not obtained. The hysteresis 
can be qualitatively explained through the coexistence region in Fig. [2j For 
U c < U < U c i the T = ground state is a FL. As one increases T for a fixed 
U, and crosses the spinodal on the right (transition line 2), a first order tran- 
sition to a paramagnetic insulating state occurs, which upon cooling does not 
transit to the paramagnetic insulating state until the left spinodal is crossed 
(transition line 1). Thus as one crosses the coexistence region and reenters via 
heating/cooling, thermal hysteresis would be obtained. The hysteretic behavior 
seen in the present theory is naturally very far removed from the rich experi- 
mentally observed hysteresis seen e.g. in V2O3. The coexistence of metallic and 
insulating islands has been experimentally observed in thin films [46], as well 
as in manganites [57], while within DMFT, where spatial inhomogeneities are 
completely ignored, the coexistence is just that of the metallic and insulating 
solutions. Experimentally, the resistivity does not increase monotonically with 
heating or cooling. Instead, multiple steps or avalanches are observed to ac- 
company the hysteretic behavior. While such details are absent in the present 
theory, nevertheless, we carry out a direct comparison of our hysteresis result 
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Figure 5: Left panel: The resistivity on a log-scale exhibits thermal hysteresis 
as a function of temperature when U c < U < U c i (see text for discussion). 
For U = 3.4i* second order transition is observed without hysteresis. For U = 
4.0i* hysteresis is not observed though the metal-insulator transition is of first 
order. Right panel: Comparison of theoretically obtained thermal hysteresis 
(U = 3.6t») with experimental observations in V2O3 [4"5] . 

with the one found experimentally in doped V2 03^45 (right panel of Fig. [5]). 
The purpose of such a comparison being to assess if the half-filled Hubbard 
model can capture at least the qualitative aspects of real materials. If it does, 
then the hope would be that a realistic theory based on finite dimensions in- 
cluding spatial inhomogeneities would be able to capture the experimentally 
observed behavior quantitatively. The best fit of the hysteresis result for a 
specific interaction (U — 3.6t«), yields t„ ~ 7305K (0.63eV). This agrees well 
with independent bandstructure calculations for V2O3 [45] and is of the right 
magnitude. 

We now focus on optical transport in the next subsection. 
3.3 Optical conductivity: 

As can be naturally expected, changes in interaction strength [//£* and espe- 
cially the Mott MIT affect dynamical or optical transport properties strongly. 
Fig. |5] shows the computed T = optical conductivity a(uj) as U/t* increases 
from a low value of 1.0 to a moderately strong value of 3.0. The inset shows 
the corresponding spectral functions. Several very interesting features can be 
seen. As U increases, a strong absorption feature emerges at ~ [7/2 for U > 2i*, 
while a second peak at ~ U emerges beyond 3i*. For low values, none of these 
features may be distinguished. The first peak arises because of excitations be- 
tween either of the Hubbard bands and the Fermi level, while the second peak 
represents excitations between the lower and upper Hubbard band. The zero fre- 
quency Drude peak is also present, but is not visible, since it has a Dirac-delta 
function form. The changes in optical conductivity are naturally understood 
through the inset of Fig. [6] For low values of U /£*, a single featureless spectral 
function is obtained, while the emergence of distinct Hubbard bands in the spec- 
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Figure 6: (Color online) The main panel shows the zero temperature optical 
conductivity for U/t* values ranging from 1.0 to 3.0 (indicated by numbers). 
The inset shows the corresponding spectral functions. 

tra mark the emergence of the first and second absorption peaks in the optical 
conductivity. At around the half bandwidth (ui ~ 1.15i*) a universal crossing 
point is is seen in the spectra. 

The effect of Mott transition on the optical absorption is illustrated clearly 
in Fig. [71 which is similar to the previous figure, except that the values of £//i* 
considered here increase from 3.0 to 5.0. The transition from the correlated 
metallic phase to the Mott insulator phase occurs at U/t* ~ 4.5, where, within 
IPT, a large gap ~ U — 2t* is known to form. The main panel shows optical 
conductivity as a function of w/t*. In the metallic phase (U < U C 2), the first 
absorption peak is seen to get narrower and surprisingly gets red-shifted as the 
Mott transition is approached. The second absorption peak begins to dominate 
as U — > U C 2 and becomes the sole feature in the Mott insulating phase. The U > 
U C 2 optical conductivity is seen to possess a clear optical gap, which increases 
with increasing U and reflects the presence of the gap in the density of states 
(see inset). 

The temperature evolution of the optical conductivities is equally interest- 
ing. In Fig. [8] we show the <t(w; T) behavior for various T's for a moderately 
strong interaction strength of U = 3.0i*. The inset again shows the correspond- 
ing spectral functions. The first absorption peak at lower temperatures loses 
spectral weight as temperature increases, which is gained by the second peak, 
and again an almost universal crossing point is seen, marking the frequency 
across which the transfer of spectral weight occurs. The Drude peak at uj = 
diminishes in height, consistent with the dc conductivity values, and disappears 
completely at T = 0.2 by forming a shoulder-like feature at higher frequency 
<jj 2.5i*. 

An earlier investigation 49 (though it was not exactly calculated for the half- 
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Figure 7: (Color online) The main panel shows the zero temperature optical 
conductivity for U/t* values ranging from 3.0 to 5.0 (indicated by numbers). 
The inset shows the corresponding spectral functions. 



filled HM) claimed that the shoulder formation is actually shifting of Drude peak 
(pseudo-Drude peak) and the phase is metallic. However, here we argue that 
this arises due to the pseudogap formation in the DoS (see the left inset) and it 
is actually an insulating state since dp(T)/dT < (see the right inset) and the 
phase lies in the crossover regime. 

The universal crossing point, known as the isosbestic point in the metallic 
phase has been observed in several materials, e.g. V2O3 [50], N^-aSe^ [51], 
La 2 - 2 ;Sr a ;Cu04 [52], Lai-^Sr^TiOa [53]. Occurrence of isosbestic points are not 
well-explained. Nevertheless such a point is believed to have a close connection 
with /-sum rules and its location is associated to microscopic energy scales in 
correlated systems [531 155] . Here, we see that the dynamics also exhibit a similar 
feature indicating an even more general basis for its existence. If we use our 
earlier estimation of t» (0.6 eV) for V2O3 or the LDA-calculated value, we find 
that the isosbestic point arises at ~ 0.9-1.0 eV which lies in the mid-infrared 
range, and which is close to that seen in a recent infra-red spectroscopy measure- 
ment [5D] i.e. at 6000 cm -1 (~ 0.7 eV). Finally we show the optical conductivity 
for the Mott insulating regime (t//i* = 4.7) evolving with temperature in Fig. [9] 
A single absorption peak is seen at low temperatures, and as T increases, spec- 
tral weight is transferred from this peak to lower frequencies, and the absorption 
peak diminishes, and experiences a blue shift. The spectral function exhibits 
negligible change as a function of temperature and although, there is indeed, 
an exponentially small rise in the density of states in the neighbourhood of the 
Fermi level, it visually appears to coincide with the T = DoS. 
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Figure 8: (Color online) The main panel shows the temperature evolution of 
optical conductivity for U/t* — 3.0. The left inset shows the corresponding 
spectral functions and the right inset shows the resistivity at T = 0.2t*. M 
marks for metallic (dp/dT > 0) and I marks for insulating (dp/dT < 0) behavior. 

4 Conclusion 

A systematic study of the half-filled Hubbard model within dynamical mean 
field theory is carried out, with the focus being on universality, scaling and 
qualitative comparison to experiments. We reformulate a well-known and ex- 
tensively employed impurity solver for the effective impurity problem that the 
lattice problem gets mapped onto within DMFT, namely the iterated pertur- 
bation theory, such that some of the problems with previous implementations 
have been overcome. We find that the coherence peak in the resistivity is a 
universal feature. A comparison with experimental measurements of resistivity 
in Se doped N1S2 with varying pressure yields qualitatively excellent agreement. 
Thermal and pressure driven hysteresis is shown to be qualitatively explicable 
within this scenario, and again, a comparison of thermal hysteresis with exper- 
iments in V2O3 are seen to yield a reasonable number for the hopping integral. 
The transfer of spectral weight across the Mott transition and the isosbestic 
points have been highlighted in the study of optical conductivity. We conclude 
that the Hubbard model does indeed represent an appropriate phenomenolog- 
ical model that can qualitatively explain a large range of phenomena observed 
in transition metal oxides. This offers hope for more detailed material specific 
studies such as those employing LDA+DMFT [56l EH [58] approaches to obtain 
quantitative agreement with experiments. 
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Figure 9: (Color online) The main panel shows optical conductivity for temper- 
atures T/t* = 0.02,0.5,0.9 and 1.5. The inset shows the corresponding spectral 
function which exhibits negligible changes as temperature is increased. 
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Appendix A: Analytical calculation of U& and TJ& 

at T = 

In this section, we demonstrate that a pole structure ansatz for the Green's 
functions and self energy combined with the IPT equations yields U c \ and U C 2, 
the bounds of the U interval where the metallic and insulating solutions coexist 
at zero temperature. 

If Q has simple poles at ±luq with equal residue a/2 (due to symmetry), then 
the singular part of Q may be expressed as 



0M = f 



i i 



UJ — CJq LU + LUq 



(18) 



Now using the Dyson equation (Eq. (|5j|) and the moment expansion of the 
Hilbert transform, we get 

2j(u;) = — uj + L(wJ 



a uj uj — 

Thus 



E(w)=w ^ (19) 

— — UJ 

a u 
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Therefore poles occur at 



= ±- 



vr 



(20) 



Using the structure of Q in the IPT equations, it is easy to see that poles at 
±Wo in Q will give rise to poles in the self energy at ±a->o and ±3u;o- Only the 
latter are physically acceptable, since the former leads to a self-consistent value 
of the residue a equal to 0. So when uj' = ±3luq, a turns out to be 8/9. 
The singular part and the FL part of the self energy may be expressed in a 
combined way as 



E(«) 



U 2 a 3 



1 



ui — 3luq 



- 3w 



This gives 



E(w ) 



2„,3 



TPa 
32wq 



1 



u> 



Now Q 1 (wo) 

E(^o) = -G- 1 ( Wo 



0. So the Dyson equation (Eq. ([5)) gives 
M 2 



-7(w ) 



and hence, after rearranging we get 



2 

w 



Z Ah 



= -u; + £(cj ) - 



32 



w - E(w ) 



(21) 



(22) 



(23) 



(24) 



The above result shows that the pole position is proportional to the square root 
of the quasiparticle weight. The transition from the metallic to the insulating 
regime occurs when ujq = O.This occurs when 
a. = 8/9, 

Uc2 



^j§£, i.e. for HCL, with 



27 



r 3 /2 4^/2 



(25) 



To get U c i, we refer back to Eq. (fT3")) which holds in the insulating phase. This 
self consistent equation for the pole residue in the insulating phase a has real 
roots only when U 2 > 27 M 2 , i.e. for HCL 



U cl 



3^3 
V2 



(26) 



The above analysis shows that we can analytically estimate U c \ and \J c i within 
the same framework. 
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